AD-A142  94?  EFFECT  OF  INERTIA  ON  FINITE  NEAR-TIP  DEFORMATION  FOR  i/1 

FAST  MGEE-1 1 1  CRCCK.  .  <U)  NORTHWESTERN  UNIV  EVANSTON  IL 
STRUCTURAL  MECHANICS  LAB  J  D  ACHENBACH  ET  AL.  FEB  84 
UNCLASSIFIED  NU-SML-TR-84-2  N00014-76-C-0062  F/G  20/11  NL 


OTIC  FILE  COPY  AD- A 142  947 


EFFECT  OF  INERTIA  ON  FINITE  NEAR-TIP  DEFORMATION 
FOR  FAST  MODE-III  CRACK  GROWTH 


J.D.  Achenbach  and  N.  Nishimura 


Department  of  Civil  Engineering 
Northwestern  University 
Evanston,  IL.  60201 


Office  of  Naval  Research 


N00014-76-C-0063 


February  1984 
NU-SML-TR-No . 84-2 


Approved  for  public  release;  distribution  unlimited 


dtic 

JUL1  1  1884Ii 


{g'4  07  10  163 


Abstract 


The  combined  effects  of  finite  deformation  and  material  inertia  have 
been  analyzed  for  fast  crack  growth  under  anti-plane  loading  conditions. 

A  steady-state  dynamic  solution  has  been  obtained  for  the  finite  strain  on 
the  crack  line,  from  the  moving  crack  tip  to  the  moving  transition  boundary 
with  the  zone  of  small  strains.  The  crack  propagates  in  a  material 
with  a  response  curve  in  uniform  shear  that  is  linear  at  small  strains, 
and  that  remains  constant  once  a  critical  strain  has  been  exceeded.  The 
corresponding  quasi-static  solution  is  given  in  the  full  zone  of  large 
deformation.  For  the  dynamic  formulation,  an  explicit  expression  for  the 
crack-line  strain  has  been  obtained  by  expanding  the  displacement  in  a 
power  series  in  the  distance  to  the  crack  line,  with  coefficients  which 
depend  on  the  distance  to  the  moving  crack  tip.  Substitution  in  the 
equation  of  motion  yields  a  nonlinear  ordinary  differential  equation  for 
the  relevant  coefficient,  which  can  be  solved  rigorously.  The  finite 
deformation  crack-line  fields  have  been  matched  to  appropriate  small-strain 
fields  at  the  transition  boundary.  The  principal  result  is  that  the 
dynamic  strain  remains  bounded  at  the  crack  tip,  apparently  due  to  the 
effect  of  material  inertia.  The  crack-line  strain  has  been  plotted  for 

f ' 

several  crack-tip  speeds.  It  decreases  with  higher  crack-tip  speed.  An 
explicit  expression  has  been  given  for  the  extent  of  the  zone  of  finite 
deformation,  as  a  function  of  the  crack  tip  speed  and  the  far-field  loading. 
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1.  Introduction 


Several  recent  investigations  have  been  concerned  with  the  effects  of 
finite  deformations  on  the  fields  of  stress  and  deformation  near  a  crack 
tip.  Most  notable  are  a  series  of  papers  by  Knowles  and  Knowles  and 
Sternberg,  in  which  full  nonlinear  equilibrium  theory  of  homogeneous  and 
isotropic,  incompressible  elastic  solids  has  been  employed  to  construct 
fields  near  the  tip  of  a  crack.  A  summary  has  been  presented  in  Ref.[l]. 

The  most  complete  results  have  been  obtained  for  a  crack  in  a  body 
subjected  to  finite  anti-plane  shear,  as  reported  in  Refs. [2]  and  [3]. 

Related  work  for  the  anti-plane  shear  case  has  been  carried  out  by  Lo  [4], 
who  found  that  for  certain  forms  of  the  strain  energy  density,  the  formu¬ 
lation  of  the  problem  becomes  formally  identical  to  those  of  certain  previously 
studied  small  strain  elastic-plastic  problems. 

The  present  paper  is  also  concerned  with  finite  deformations  near  a 
crack  tip  under  anti-plane  shear  loading.  The  formulation  of  the 
governing  equations  follows  by-and-large  Refs. [2]  and  [3].  A  notable 
difference  is,  however,  that  we  study  a  growing  crack,  and  that  the  effects 
of  inertia  have  been  included  in  the  formulation.  We  consider  a  strain 
energy  density  which  yields  a  response  curve  for  uniform  shear  that  is 
linear  at  small  strains,  but  that  maintains  a  constant  plateau  once  a 
certain  critical  strain  has  been  exceeded. 

The  method  of  solution  employs  power  series  expansions  in  the  distance 
to  the  crack  line,  with  coefficients  which  depend  on  the  distance  to  the 
crack  tip.  To  analyze  small  strain  crack  line  fields  in  elastic  perfectly- 
plastic  materials,  such  expansions  have  previously  been  used  by  Achenbach 


2 


and  Dunayevsky  [5]  and  Achenbach  and  Li  [6]  for  quasi-static  problems, 
and  by  Achenbach  and  Li  £ 7 )  for  dynamic  crack  growth  cases.  Substitution 
of  a  displacement  expansion  in  the  equation  of  motion  yields  a  nonlinear 
ordinary  differential  equation  for  the  relevant  coefficient,  which  can  be 
solved  rigorously.  The  crack-line  fields  have  been  matched  to  appropriate 
small  strain  elastic  fields  at  a  transition  boundary. 

The  paper  also  includes  a  section  which  shows  that  for  dynamic  crack 
growth  the  stresses  and  the  displacement  gradients  are  continuous  across 
the  moving  transition  boundary  which  precedes  the  crack  tip,  at  least  near 
the  crack  line.  Continuity  of  these  quantitites  all  around  the  transition 
boundary  has  been  shown  for  the  quasi-static  case.  Furthermore,  the  full 
quasi-static  field  has  been  obtained  in  the  zone  of  large  deformation.  For 
the  elastic  material  that  is  considered  in  this  paper,  the  static  fields  for 
a  moving  and  a  stationary  crack  are  identical  for  coordinates  centered  at  the 
crack  tip.  The  quasi-static  crack-tip  strain  is  shown  to  be  singular.  An 
Independent  application  of  the  expansion  method  yields  the  correct  crack-line 
approximation  to  the  exact  solution. 

The  analysis  of  the  crack-line  fields  for  steady-state  dynamic  crack 
growth  has  yielded  the  most  surprising  result  of  this  paper.  It  is  shown 
that  the  effect  of  inertia  is  to  remove  the  strain  singularity  at  the  crack 
tip.  Thus,  for  the  dynamic  problem  both  the  stress  and  the  strain  are 
bounded  at  the  crack  tip.  Curves  are  presented  which  plot  the  strain  versus 
distance  along  the  crack  line.  An  expression  has  been  obtained  for  the 
length  of  the  zone  of  large  deformation  as  a  function  of  the  crack-tip  speed. 
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2 .  Governing  Equations 

The  equations  which  govern  finite  anti-plane  shear  deformation  in  an 
unbounded  body  containing  a  crack  have  been  formulated  by  Knowles  f  2 1  and 
Knowles  and  Sternberg  [3].  In  this  Section  we  briefly  review  the  formulation 
of  [2]  and  f3]  in  a  slightly  different  notation,  we  introduce  a  different 
strain-energy  density  than  was  used  in  these  papers,  and  we  extend  the 
formulation  to  the  dynamic  case. 

Let  and  x^  be  the  Cartesian  coordinates  of  a  typical  point  in  the 
undeformed  and  deformed  states,  respectively.  For  anti-plane  strain  we 


x_  -  X,  +  U(X, ,X„,t),  x  -  X 
3  3  12  a  ct 


a  =  1,2  . 


(2.1) 


The  corresponding  components  of  the  deformation  gradient  F  are 


(2.2) 


and  the  first  fundamental  scalar  invariant  is 


Ix  -  tr(F  F)  -  3  +  e2. 


where 


(2.3) 


e  *  U  U  ,  U  -  3U/3X 
,a  ,a  ,a  a 


(2.4a,b) 


We  consider  an  incompressible  material.  Hence  the  deformation  is  locally 
volume  preserving,  and  the  Jacobian  of  the  mapping  (2.1)  is 


J  -  det(F)  -  1 


(2.5) 


Let  t  be  the  actual  (Cauchy)  stress  tensor  field  in  the  deformed 


body,  and  let  S(Xa>t)  be  the  corresponding  nominal  (Piola)  stress  field 

in  the  undeformed  configuration.  Then 

T  T  —1 

t  -  S  F  ,  S  -  t(F; 


(2.6a,b) 


The  analysis  of  this  paper  is  restricted  to  the  class  of  incompress¬ 


ible  solids  for  which  W  depends  on  1^  only 
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W  -  Wd^  for  all  I  _>  3,  W(3)  =  0. 


As  discussed  by  Knowles  [1],  we  then  have 


t3o  *  '.3  S3o  Sa3 


2"'<I1>U,„- 


V  *  Sa6  *  S33  '  0  •  t33  ’  2W'(Il>e2 

where  e  is  defined  by  (2.4a,b),  and 
W’d^  -  dW/dl^ 

For  dynamic  problems  these  equations  are  supplemented  by 
S3a,o  "  poS  •  * 

where  pQ  Is  the  mass  density  In  the  reference  configuration. 

To  define  the  desired  relation  between  S_  and  U  ,  we  consider  the 

3a  ,a 

static  problem  of  an  unbounded  body,  which  Is  subjected  to  the  following 

condition  at  Infinity 

U(X.  ,X_)  -  kX_  +  od)  as  XX  -  »  . 

12  2  a  a 


(2.7) 

(2.8) 

(2.9) 

(2.10) 

(2.11) 


The  solution  then  is 


U(X1,X2)  -  <X2 


for  all  (XX,X2) 


T23  "  T32  "  *  2W’(3+k2)k,  t33  -  t(k)k  (0  <_  k  <  ») 

The  graph  of  t  versus  k(0  £  <  <  •)  is  called  the  response  curve  in 
simple  shear.  In  this  paper  we  wish  to  consider  a  response  curve  of  the 
form  shown  in  Fig.  1,  which  is  defined  by 


(2.12) 


(2.13) 


(2.14a  ,b) 


TOO  - 


UK  K  <  E 
—  O 


ue  k  >  e 
o  —  o 


(2.15a) 


(2.15b) 


The  corresponding  strain  energy  density  for  an  incompressible  elastic 
solid  then  follows  from  (2.14a)  as 


[i>  -  / 


d1-3> 


t (<)d<  (3  _<  1^  <  °°) 


Substitution  of  (2.15a,b)  into  (2.16)  yields 


w(i1) 


f*«r3> 

.“KV3)eollS  - 


3  <  I,  <  3  +  ez 
—  1  —  o 


3  +  e  <  I 


o-l 


The  stress-strain  relations  follow  from  (2.8)  and  (2.17)  as 


S..  «  wU  „ 
3a  ,a 


for  0  <  c  <  e  • 
—  —  o 


S  *  S  U  /e  for  e  >  e  , 
3a  o  ,a  —  o 


where 

S  *  ye 
o  o 

It  is  noted  that  (2.19)  implies  that  for  e  >  we  have 

sh  +  sh’sl  ■ 

It  is  noted  that  the  Cauchy  stresses  do  not  satisfy  a  relation  of 


analogous  simplicity. 


(2.16) 


(2.17a) 


(2.17b) 


(2.18) 


(2.19) 


(2.20) 


(2.21) 


4.'  t.'k'  O 


ijL  tttttTjTTTv  ’*  '■'  77?7V  7  T 


3.  Discontinuities  across  a  Propagating  Transition  Boundary 
A  surface  which  separates  regions  where  the  constitutive  behavior  is 
defined  by  (2.18)  and  (2.19),  respectively,  will  be  called  a  transition 
boundary.  In  this  Section  we  examine  the  possibility  of  propagating 
discontinuities  in  the  stress,  S^(X,t),  the  displacement  gradient, 

Fjj  (X, t) ,  and  the  particle  velocity,  x^X.t),  across  a  propagating 
transition  boundary,  for  the  case  that  the  linear  relation 

(2.18)  holds  ahead  of  the  moving  surface,  while  the  nonlinear  relation 

(2.19)  applies  behind  it. 


The  propagating  surface  is  represented  by  the  relation 

Z(X,t)  -  0  (3.1) 

The  unit  normal  to  the  surface  is  denoted  by  N.  The  speed  of  propagation 
VN»  which  is  the  speed  with  which  the  surface  traverses  the  material,  is 
positive  when  the  surface  moves  in  the  direction  of  A  discontinuity  of  a 
field  quantity,  say  G(X,t),  is  denoted  in  the  usual  manner  by 

[ [G] ]  =  G+  -  G"  (3.2) 

where  G+  and  G  are  the  limits  as  X  approaches  a  point  on  the  surface  along 
paths  entirely  ahead  of  and  entirely  behind  it  . 

The  displacement  is  continuous,  and  hence  displacement  gradients  in 
the  surface  are  also  continuous 


[[3xi/arg]]  -  0 


f 


B  -  1,2 


(3.3) 


where  and  T ^  are  two  independent  coordinates  in  the  surface  I(X,t)  =  0.  The 
balance  of  linear  momentum  yields  the  well-known  relation 

t[SiJ13Nj  *  ~poV[ii]1  *  (3.4) 

where  is  the  mass  density  in  the  reference  configuration,  and  x 
denotes  the  particle  velocity.  From  the  balance  of  energy  it  follows 
that  the  energy  dissipated  locally  as  the  surface  passes,  may  be  written 
as 

H  -  [  [W  +  |  P£)xiii]  ]VN  +  [  [xiSij  ]  ]N  ,  (3  5) 

where  W(X,t)  is  the  strain  energy  density.  It  is  assumed  that 

H  >_  0  ,  (3.6) 

i.e.,  the  energy  dissipation  across  a  propagating  transition  boundary 
cannot  be  negative  [  8  ] . 

It  can  be  verified  that  H  may  be  rewritten  as 


I 


we  finally  find 


H  =  {  [  [W]  ]  -  [[3xi/3N]]  (S^  +  S^)^}  VN 

For  anti-plane  shear,  (3.11)  reduces  to 

H  -  {[[WH  -  |  [[3x3/3N]](S^  +  S3a)Nc(}VN 
where  a  **  1,2.  An  expression  similar  to  (3.12)  has  been  presented  in  [9] 


In  this  paper  the  interest  is  focussed  on  a  field  x^X^.X^t)  which  is 
antisymmetric  with  respect  to  the  X^  axis,  but  which  has  discontinuities 
across  a  transition  boundary  which  is  symmetric  relative  to  X^  =  0,  and  on  which 
>  0.  This  case  applies,  for  example,  to  a  transition  boundary  ahead  of  a 
Mode-Ill  crack  tip.  The  condition  H  ^  0  then  becomes 


'H 


H  -  H/Vn  -  [ [W] ]  -  |[[3x3/3N]](S3a  +  S3a)NQ  >  0 

Let  us  first  consider  the  intersection  point  on  X2  *  0,  which  is  defined 
by  ■  0,  X^  ■  XQ(t).  Because  of  antisymmetry  we  have  x3(X^,0,t)  =  0, 
and  hence 

3x3/3X1  h  0  , 

while  (3.3)  implies 

e  -  e+  ■  3x-/3X-  “  e 

J  2  O 

where  e  is  defined  by  (2.20).  By  using  the  expressions  for  W  given  by 
(2.17a)  and  (2.17b),  it  then  follows  that 

H  -  t  [W]  ]  -  -|v(e+)  -  ye~eo  +  -  0 


(3.13) 


3. 


Next  we  consider  a  position  (X^dX^) .  Here  we  have  by  (3.13)  and 
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dH  >  0  , 


where 


dH  *  uc+de+  -  yeod e"  -  \u+  ~  e")  (yde*  +  — °  deN 

e 


ye  e 
°  N 

(O2 


de  ) 


1  ,.  +  ,  -w  +  Nv 

*2  <dEN  '  deS)<UEll  +  U  o  ~=> 

£ 


In  (3.18), 

eN  -  3x3/3N 

Since  e+  *  e  ■  e  and  e*  »  e_ *  0  ,  (3.18)  reduces  to 

O  N  N 

peo(d£+  -  de  )  _>  0 


(3.17) 


(3.18) 


(3.19) 


(3. 2D) 


In  the  domain  of  small  deformation  the  quantity  e  cannot  increase  with  distance  from 

the  crack  line,  hence  de+  <_  Q.  On  the  other  hand  e  maintains  the  condition  of 
large  deformation,  e  _>  £q,  behind  the  transition  boundary,  and  hence  de  ^.0. 
Consequently  (3.20)  can  be  satisfied  only  if 

de+  -  de".  =  0  (3.21) 

Equations  (3.15)  and  (3.21)  imply 

[[e]]  -  0  (3.22) 

Equation  (3.22),  together  with  the  constitutive  equations,  (2.18)  and  (2.19), 
give 


[[SjJ]  -  y[[3x3/3Xo]]  -  wNa[[3x3/3N]]  (3.23) 

where  (3.3)  has  also  been  used.  Next  (3.4),  (3.10)  yield 

ttS3allNa  "  P0VNt[3x3/3N3]  “  0  (3.24) 

Substitution  of  (3.23)  subsequently  gives 


U(1  -  Mj)  [[3x3/3N]]  *  0 


where  M*  =  p  Vz/u  .  It  follows  that 
N  ON 

[[3x3/3N]]  =  0 


(3.25) 


(3.26) 


provided  that  M^<  1.  Equation  (3.26)  in  turn  implies  by  (3.10)  that 
=  0  ,  and  hence 

t[S3a]J  =  0  (3.27) 

by  (3.23).  We  have  thus  proven  the  continuity  of  the  relevant  quantities 

across  the  transition  boundary  near  the  crack  line,  for  dynamic  problems. 

In  (quasi)  static  problems  (p  =  0) ,  we  can  prove  the  same  results  all  along 

o 

the  transition  boundary.  Indeed,  (3.3)  and  (3.4)  give  [10] 


[[S3a]][[3x3/3Xa]]  =  0  , 
which  is  written  explicitly  as 

I  =  {y(9x,/9X  )+  -  S  (9x_/9X  )~/e~}{ (9x_/9X  )+  -  (9x,/9X  )"}  =  0 

3  a  o  3  a  3a  3a 

It  follows  that 

+  4-  —  —  +  —  — 

I  *  y(e  )  +  S  e  -  S  (9x,/9X  )  (9x,/9X  )  /e  -  u(9x,/9X  )  (9x,/9X  )  . 

o  o  3  a  3a  3a  3a 


By  using  the  inequality 

,  +  -  .  ,  +  +  -  -  sk 

lx3,ax3>J  -  (x3,ax3,«)  (x3,Sx3.6^  ’ 

it  follows  from  (3.30)  that 


2  2 

I  >  y(e+)  +  S  e  -  S  e+  -  ye+e  ■  u[(e+)  +  e  e  -  e  e+  -  e+e  ] 

—  0  0  o  o 

-  li(e+  -  e”)(e+  -  e  )  >  0  , 

o  — 


(3.28 ) 


(3.29) 


(3.30) 


(3.31) 


(3.32) 


since  e  <  e  and  e  <  £  But  (3.28)  and  (3.32)  will  be  in  conflict, 

—  —  o 

unless  e+  =  e  on  the  transition  boundary  (note  that  e+  =  e  will  lead  to 
o 

e+  =  eQ  also).  It  can  be  shown  that  this  result,  together  with  (3.3)  and 

(3.4),  yield  either  e+  *  t  ,  or  N  *  Vx^  *  0  (i.e.  x^  is  locally  constant) 

on  the  transition  boundary.  As  will  be  seen  later,  we  may  exclude  the 

latter  possibility  for  the  present  application.  We  thus  conclude  that 

e+  =  e  =  e  on  the  transition  boundary.  With  this  result,  we  can  follow 
o 

the  procedures  in  Eqs.(3.23)  -  (3.27)  to  prove  continuity  of  S^.Sx^/aN  and 


Static  Fields  Near  a  Crack  Ti; 


In  the  immediate  vicinity  of  a  crack  tip  the  deformation  is  large, 
i.e.,  e  >  cq  ,  and  consequently  Eqs.(2.19)  and  (2.21)  hold.  Outside  the 
region  of  large  deformation,  linearized  elasticity  according  to  Eq.(2.18) 
applies.  The  two  regions  are  separated  by  a  transition  boundary.  In  this 
section  it  is  assumed  that  the  largest  length  dimension  of  the  region  of 
large  deformation  is  small  as  compared  to  the  length  of  the  crack.  This 
assumption  implies  that  the  crack  may  be  taken  as  semi-infinite,  and  that 
the  displacement  outside  the  region  of  large  deformation  may  be  represented 
by  the  well-known  square-root  solution  according  to  linearized  elasticity. 

The  conditions  on  the  faces  of  the  crack  are 

U  2(X1*0±)  “  0  —  <  Xx  <  0 

The  equation  which  governs  the  displacement  field  in  the  zone  of  large 
deformation  is  obtained  by  substituting  (2.19)  into  the  balance  equation 
for  static  equilibrium.  The  result  is 


(U 


2U,l",2°,21  +  (U,1>2U,22 


0  . 


Equation  (4.2)  can  be  rewritten  as 


(1/U  .)L{U  .}  -  (1/U  ,)L{U  -  0  , 

>1  *1 

where  the  differential  operator  L  is  defined  as 


L  -  (3/3X,)  -  (U  ./U  ,) (3/3X_) 


It  is  now  noted  that  L  is  a  total  differential  along  the  characteristic 
line  defined  by 

dX2/dX1  =  -  Ufl/U>2  (4.5) 

Along  such  a  line  (4.3)  can  be  integrated  to  yield 

£n(U  ./U  .)  -  C  ,  or  U  ./D  «  -D  ,  (4.6a,b) 

where  C  and  D  are  constants.  By  combining  (4.5)  and  (4.6b),  we  conclude 
that  the  characteristics  are  straight  lines: 

dX2/dXx  -  D  (4.7) 

On  these  characteristics  the  displacement  U  is  constant,  because  it  follows 
from  (4.5)  that 

dU  -  0  (4.8) 

The  problem  formulation  in  the  zone  of  large  deformation  as  defined  by 
Eqs.(2.21)  and  (4.8)  is  now  identical  to  the  one  for  small  strain  deformation 
near  a  crack  in  an  elastic  perfectly-plastic  material  with  the  Huber-Mises 
yield  condition.  The  solution  to  the  latter  problem  was  given  first  by  Hult 
and  McClintock  [11].  It  was  discussed  in  some  detail  by  Rice  [12].  In  the 
present  context  the  Hult-McClintock  solution  shows  that  the  zone  of  large 
deformation  is  a  circle  of  radius 

Ro  ■  h  <Wso>2  <4-9> 

which  is  centered  at  a  point  E  on  the  crack  line  (see  Fig.  2),  a  distance 


R 


o 


ahead  of  Che  crack  tip. 
/R  2 

U  *  {  2?j  p  *111  Sin9 


In  the  zone  of  large  deformation  we  have  [11], [12]: 

(4.10) 


*  -S  sin8  ,  S.0  ■  S  cos0  (4.11a,b) 

jl  O  51  O 

In  the  region  of  linearized  elasticity  we  have 

u  ■  (ht  u  hn  sl4» 

h 

531  ■  •  fe)  *111  sl4* 

532  “  ( 2itr)  \ll  ““T*1 

where  (R,i]»)  are  polar  coordinates  centered  at  the  point  E.  Thus,  (4.10)  - 

(4.11a,b)  hold  for  R  <  R  ,  while  (4.12)  -  (4.14)  are  valid  for  R  >  R  . 

o  o 

We  conclude  this  section  by  rederiving  the  field  on  the  crack  line  in 
the  zone  of  large  deformation,  by  the  use  of  an  expansion  method.  This 
method  can  be  used  for  problems  that  are  not  otherwise  solvable.  The 
method  is  shown  here  to  demonstrate  its  utility,  and  as  a  preliminary  to 
its  application  in  the  next  section  to  the  problem  of  dynamic  crack  growth. 

Following  the  approach  of  Achenbach  and  Dunayevsky  [5]  and  Achenbach  and 
Li  [6],  we  seek  a  solution  to  Eq.(4.2)  for  small  values  of  X2(X2/X^  <<  1), 
in  the  form 


(4.12) 

(4.13) 

(4.14) 


U  -  U(1)(X,)X„  +  U(3)(X,)X^  +  0(X*) 


(4.15) 
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Substitution  of  (4.15)  into  (4.2)  followed  by  collection  of  terms  of 
order  yields 

u(l)u(1)  _  2(U(J’))2  -  0  (4.16) 

Equation  (4.16)  is  satisfied  by 

-  (1/C)(U(1))2  (4.17) 

where  C  is  a  constant.  The  general  solution  to  (4.17)  is  =  -C/(X^+D). 

A  singular  and  multiple-valued  solution  at  the  crack  tip  is  obtained  by 
setting  D  0.  Hence 

U(1)  -  C/Xj^  (4.18) 

where  C  is  a  still  to  be  determined  constant.  At  X.  ■  X  ,  X_  ■  0,  we  have 

1  o  z 

■  e  ,  and  thus  C  -  X  e  where  X  is  the  X,  coordinate  of  the  transition 
o  o  o  o  1 

boundary  on  the  crack  line.  It  follows  that 

U  -  (X /X  )e  X  +  O(x’)  (4.19) 

o  i  o  i  i 

The  corresponding  stresses  at  the  transition  boundary  follow  from  (2.19) 

as 

S31  -  -  So(X2/Xo)  ,  S32  -  So  (4.20a,b) 

Near  the  crack  line  in  the  zone  of  linearized  elasticity,  we  now  assume 
solutions  of  the  general  form  (4.12)  -  (4.14).  Continuity  of  S31  and  S32 
then  yields  by  the  use  of  (4.13)  and  (4.14) 

9  *  X2/Xq  -  I*  (4.21) 


.  r.  r. 


rvv.Tw 


which  implies  that  X  ■  2R  .  It  is  now  easily  verified  that  for  6  <<  1, 

o  o 

(4.19)  and  (4.20a,b)  do  indeed  correspond  to  (4.10),  (4.11a,b)  with  (4.9), 
respectively. 

As  a  final  comment  we  note  that  the  results  derived  in  this  section 
also  apply  to  the  crack  line  fields  for  quasi-static  crack  growth.  For  a 
propagating  crack  it  is  convenient  to  define  the  coordinate  system  (X^.X^) 
as  moving  with  the  crack  tip.  Stationary  coordinates  (X,Y)  and  moving 
coordinates  (X^.X^)  are  related  by 

^  -  X  -  a(t)  ,  X2  *  Y  ,  (4.22a,b) 

where  X  ■  a(t),  Y  *  0  defines  the  position  of  the  crack  tip  in  the  (X,Y) 
system.  In  the  moving  (X^,X2)  system,  the  material  time  derivative  is 
defined  as 

C)  -  (3/3t)  -  a(3/3X1)  ,  (4.23) 

where  a  -  da/dt  is  the  crack-tip  speed.  Spatial  derivatives  are  the  same  in 
the  two  systems.  Since  the  time  derivatives  do  not  appear  in  the  elastic 
constitutive  equations  used  in  this  paper,  they  do  not  enter  the  formulation 
of  the  quasi-static  crack  growth  problem.  Hence,  the  governing  equations 
and  (therefore , the  fields  are  identical  for  the  stationary  and  the  moving 
crack  tip,  provided  that  they  are  defined  relative  to  a  coordinate  system 
that  is  centered  at  the  (stationary  or  moving)  crack  tip. 
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5.  Dynamic  Fields  near  a  Moving  Crack  Tip 
It  is  convenient  to  formulate  the  governing  equations  in  the  moving 
coordinates  (X^,X£)  defined  by  (4.22a,b),  with  the  material  time  derivative 
as  defined  by  (4.23).  In  this  paper  we  will  restrict  the  attention  to  the 
steady-state  field  near  a  crack  tip  which  moves  with  a  constant  speed  c. 

The  material  time  derivatives  then  reduce  to 

(*)  -  -cO/dXj),  (  )  -  c2(32/3X2)  (5 .la,b) 

Relative  to  the  moving  coordinate  system,  the  equation  of  motion  (2.11) 
becomes 


S-  r.  =  P„c  u  n  •  where  a  -  1,2 
3a, a  o  ,11 

Substitution  of  (2.19)  into  (5.2)  yields  in  the  zone  of  large  deformation 


(5.2) 


-  2B,1B,2B,12  +  <U,1)!U.22  *  X[ <B,l)a+(U.2)J,3/2u.n 


(5.3) 


where 

X  -  c2p  /S  (5.4) 

o  o 

A  solution  in  the  vicinity  of  the  crack  line  can  be  obtained  by  using 
the  expansion  given  by  Eq.(4.15).  By  substituting  (4.15)  into  (5.3)  and 
collecting  terms  of  order  X^  we  find 

u(l)u(1)  _  2(U(J))2  -  X(U(1))VjJ  (5.5) 

For  X  ■  0,  (5.5)  reduces  to  (4.16).  It  is  convenient  to  define 


(5.6) 


and  to  rewrite  (5.5)  in  the  form 


dU(1)  _  dV^ 

U(1)(1-XU(1))  2V(1) 

This  equation  can  be  integrated  to  yield 


(5.7) 


(5.8) 


where  Vq  is  a  positive  constant*  and  the  sign  is  chosen  so  as  to  make  the  right 
hand  side  positive.  For  a  few  values  of  X  *  ±V^^/V  has  been  plotted  in 
Fig.  3.  Also  plotted  in  Fig.  3  is  the  corresponding  relation  for  X  =  0, 
which  follows  from  (A .19)  and  (5.6)  as 


e  X  V(1) 
o  o  qs 


-  ^ 
qs 


(5.9) 


Since  we  expect  that  U 


(1) 


U(1>  and  V(1> 

qs 


as  X  -*■  0  ,  it  can  be 
qs 


seen  from  Fig.  3  that  the  relevant  dynamic  solution  for  and  should 

be  located  in  the  part  of  the  phase  plane  defined  by 


0  <  U(1)  <  1/X  ,  V(1)  <  0 

Hence,  for  Vq  >  0  ,  the  minus  sign  applies  in  Eq.(5.8).  An  Important 
consequence  of  (5.10a)  is  that  remains  bounded  for  X  +  0. 

Ue  now  return  to  Eq.(5.6),  and  we  write 

dXl  1  (U(1)-  1/X)2 

V  (U^) 
o  ' 

where  (5.8)  has  also  been  used.  Integration  yields 


(5.10a,b) 


(5.11) 


X2  U 


1 

(1) 


+  C 


(5.12) 


This  expression  contains  2  constants,  V  and  C,  which  can  be  obtained  from 

o 

the  conditions  at  X.  =0  and  X,  =  X  ,  where  X  is  the  intersection  of  the 

1  1  o  o 

transition  zone  with  the  X^  axis  in  the  moving  coordinate  system.  We  have 

U(1)  -  e  at  X  -  X  (5. 

o  1  o 

The  only  place  on  the  crack  line  where  singularities  can  be  expected  is 
X^  *  0.  Since  Eq.(5.8)  shows  that  -*■  «  as  ■+•  1/A,  we  conclude 


U(1)  =  1/A  at  X.  -  0  (5. 

Application  of  (5.13)  and  (5.14)  to  (5.12)  yields 

C  -  (2/AV  )2nA  (5. 

o 

V  -  (-M2  +  22. nM2  +  1/M2) /XX  (5.! 

o  o 

where  M  is  the  Mach  number 

M2  -  c2/(p/pq)  (5.: 


Substitution  of  (5.15)  and  (5.16)  into  (5.12)  yields 

fl  M  1-M4U2  +  2M2U&n(M2U)  (5  . 

Xo  (1-M4  +  2M22nM2)U 

Here  U  is  the  normalized  strain 

U  -  U(1)/e  (5.] 

o 

The  condition  that  the  crack  tip  is  in  the  zone  of  large  deformation 
implies  1/X  >  eQ,  and  hence  M2  <  1,  since  M2 


*  e  X  . 
o 
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V. 


ii 

•  Vj 

■V* 

i 


Plots  of  U  versus  X, /X  can  be  obtained  from  (5.18).  Results  are 

1  o 

shown  in  Fig. 4,  for  various  values  of  M.  For  M  >  0,  U  remains  finite 
(and  equal  to  1/M2)  at  X^  ■  0.  In  the  limit  M  -*•  0  we  find 

X./X  »  1/U  ,  or  U(1)  *  (X  /X. )e  ,  (5.20) 

1  o  o  1  o 

in  agreement  with  the  quasi-static  results  given  by  (4.19). 

Equation  (5.18)  is  the  main  result  of  this  paper.  It  is,  however, 

desirable  to  determine  an  explicit  expression  for  Xq  in  terms  of  the  crack 

tip  speed  and  the  far-field  loading,  albeit  by  approximation.  In  Section  3 

it  was  shown  that  3x_/3N  (or  3U/3N),  and  S,  are  continuous  at  the  transition 

3  3a 

boundary.  Thus,  the  fields  of  large  deformation  can,  in  principle,  be  matched 
to  corresponding  linearly  elastic  fields,  to  yield  an  expression  for  Xq. 

Under  the  assumption  that  the  zone  of  large  deformation  is  small,  the  small 
strain  fields  are  taken  as  the  asymptotic  near-tip  fields  for  a  crack  in  a 
linearly  elastic  solid.  For  the  quasi-static  case  the  solution  (4. 12) -(4. 14) 
shows  that  exact  matching  can  be  achieved,  provided  that  the  center  of  the 
small  strain  field  is  shifted  to  a  point  on  the  crack  line  ahead  of  the 
actual  crack  tip.  For  the  dynamic  problem  we  do  not  expect  that  the  zone  of 
large  deformation  will  be  circular,  and  hence  it  is  unlikely  that  a  simple 
solution  valid  along  the  complete  transition  boundary  can  be  obtained.  In 
any  event  we  have  expressions  for  the  large  deformation  fields  only  in  the 
Immediate  vicinity  of  the  crack  line.  It  is,  however,  noted  that  to  first 
order  in  X^  the  small-strain  quasi-static  and  steady-state  dynamic  solutions 
are  both  of  the  general  form,  given  by  (4.12)  -  (4.14).  Hence  it  seems 


reasonable  to  assume  linearly  elastic  fields  of  the  type  (4.12)  -  (4.14), 
at  least  near  the  crack  line,  but  where  should  now  be  interpreted  as 

an  elastodynamic  Mode-Ill  stress  intensity  factor.  For  a  semi-infinite 
crack  the  elastodynamic  and  quasi-static  stress-intensity  factor  are 
related  by  "  (1— M)^  (K^.^)^,  see  e*8*  (13, P.35].  It  follows  that 

(4.9)  also  holds  for  the  dynamic  problem,  but  with  the  understanding  that 
Rq  is  the  local  radius  of  curvature  of  the  transition  boundary  at  its 
intersection  with  the  crack  line. 

Near  the  crack  line  continuity  of  311/ 3N  implies  continuity  of  V^. 


At  the  transition  boundary  we  then  have 


(5.21) 
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